Long-term classical and general relativistic effects on the radial 
velocities of the stars orbiting Sgr A* 



Lorenzo lorio 1 

Ministero dell'Istruzione, dell'Universita e della Ricerca (M.I.U.R.) 
Fellow of the Royal Astronomical Society (F.R.A.S.). 
Viale Unita di Italia 68, 70125, Bari (BA), Italy 

lorenzo . iorio@libero . it 

Received ; accepted 



-2 - 



ABSTRACT 



We analytically work out the cumulative, i.e. averaged over one orbital rev- 
olution, time variations (v p ) of the radial velocity v p of a typical S star orbit- 
ing the Supermassive (M. rs 1O 6 M ) Black Hole (SBH) hosted by the Galactic 
Center (GC) in Sgr A*caused by several dynamical effects. They are the gen- 
eral relativistic gravitoelectromagnetic (GEM) fields of the SBH, its quadrupole 
mass moment Q2, and a diffuse dark matter distribution around the SBH. All of 
them induce non-zero long-term radial accelerations proportional to the eccen- 
tricity e of the orbit. By taking the S2 star, orbiting the SBH along a highly 
eccentric (e = 0.8831) ellipse with a period Pb = 15.9 yr and semi-major axis 
a = 1031.69 au, we numerically compute the magnitudes of its radial acceler- 
ations. The largest effects are due to the general relativistic Schwarzschild-like 
gravitoelectric (GE) field, with /^ GE) \ = 8.2 x 10~ 5 m s" 2 , and the diffuse 
material distribution, modeled with a Plummer-type mass density profile, with 
(v p dm ^ = 3.8 x 10~ 6 m s -2 . The effects caused by the general relativistic Kerr- 

type gravitomagnetic (GM) field and by Q2 are smaller by orders of magnitude. 
By assuming an uncertainty in measuring the radial velocities of about 15 km 
s _1 , the future accuracy in measuring (v p ) can be evaluated to be of the order 
of 2.4 x 10 -5 m s™ 2 over an observational time span At = 20 yr. Currently, the 
available radial velocity measurements cover just 7 yr. 

Subject headings: black hole physics-Galaxyxenter-relativity-techniques: radial 
velocities 
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Introduction 



It is no w widely accepted that the Galactic Cen t er (GC) hosts a S upermassive Black 



Hole (SBH) ( IWollman et all 119771 : iGenzel et al.lll996t iGhez et al.ll2008l : ISch5del et all 12002 
MeliaJ 20071) . whose p osition coincides with that of the radio-source Sagittariu s A* (Sgr 



A*) flReid et all 120071). The SBH has a mass of the order of M. = 4 x 1O 6 M (IGhez et al 



20081 ; iGillessen et all l2009al Jbl) and, consequently, a Schwarzs child radius 



0.084 au. 



In its immediate vicinity a number of rapidly orbiting stars ( iPaumard et al.l 120061 ) have 
been detected and tracked in the infrared since 1992 at the 8.2-m Very Large Telescope 
(VLT) on C erro Par anal, Chile and the 3.58-m New Technology Telescope (NTT) on La 
Silla, Chile (lEckart &: Genzellll996l) . and since 1995 at the Keck 10-m telescope on Mauna 
Kea, Hawaii (IGhez et al.lll998l ). They are mostly main-sequence stars of spectral class 
B, and are naturally used as test particles for the gravitational potential in which they 
move. Particularly important is the bright star S2 (SO- 2 in the Keck nomenclature), of 
spectral type BO-2 V, orbiting the S BH in 15.9 yr along an orbit with ellipticity e = 0.8831 
and semimajor axis a = 1031.69 au (IGillessen et al.ll2009a| ). Indeed, a complete o rbital 
revolution of it is now covered by the available data records ( Gillessen et al.ll2009al Jbl). 



If we look at the smallness of the ratio r of the average distance (r) = a(l + e 2 /2) to the 
Schwarzschild radius T g CIS Clll index of the importance of the Einstein's General Theory of 
Relativity (GTR) in several astronomical and astrophysical systems, it can be easily realized 
that the Galactic SBH and its stars is, in principle, an ideal local laboratory to put on 
the test GTR and other alt ernative theories of gr avity. Indeed, by considering the Earth's 
artificial satellite L AGEOS (ICohen Sz Smithlll985f) . the Sun's planet Mercury, the extrasolar 
planet WASP-19b (THebb et al I l2010h. the double binary pulsar PSR J0737-3039A/B 
( IBurgay et al.ll2003f ) and the star S2 in Sgr A*, we have 



^LAGEOS 



1.383319 x 10 



t'Mercury 
tWASP-19b 



2.0023 x 10 7 



8.74 x 10 5 



(1) 



tPSR J0737-3039A/B 



1.15 X 10 5 , 



tS2 



1.7 x 10 4 . 



It can be noted that r for S2 is one order of m agnitude smaller than that for PSR 



J0737-3039A/B. In this respect, several authors ( 


Jaroszvnski 


1998; 


Fragile & Mathews 


2000; 


Rubilar & Eckart 


2001; 


Weinberg et al 


. 2005: 


Kraniotis 


20071; 


Nucita et al. 


2007; 


Will 


2008; 


Preto & Saha 


2009; 


Kannan & Saha 


20091: Merrit et al 


. 2010) worked out with a 



variety of techniques and approximations the direct effects of GTR on different quantities 
characterizing the orbital motions about the SBH, mainly some Keplerian orbital elements. 



-4- 



If, on the one hand, they are useful to give "intuitive" valuable insights about the magnitude 
of the relativistic effects occurring in such a scenario, on the other hand they are not 
directly measurable. 

Concerning the motion of the S stars around the SBH in the GC, the directly observable 
quantities are the astrometric measurements of their positions in the sky in terms of right 
ascension a and declinatio n 5, and their rad i al vel ocities v p . Concerning the first kind of 
observations, according to lEisenhauer et al.l (120091 ). future astrometric measurements of 
S2 may bring its relativistic perinigr iconic precession (u m ) into the measurability domain. 
Indeed, the perinigricon rate would be indirectly inferred from the corresponding apparent 
position shift. To this aim, it must also be considered that such a shift is not as easily 
detected as it may seem since it needs to be measured from t he same data from which also 
the orbital elements have to be determined (IFritz et al.l 120101 ) . Thus, here we will focus on 
the dynamical effects directly caused by GTR and other competing classical forces on the 
radial velocity v p . On the one hand, it will be possible to straightforwardly work out in 
an analytical way the net time variations of it averaged over one orbital revolution. This 
allows for a more direct and unambiguous confrontation of the theoretical predictions with 
the observations. On the other hand, from a practical point of view the radial velocity data 
are easier to handle with respect to the astrometric observations. Indeed, the inclusion of 
new data into pre-existent records nee ds no special car e because the radial velocities refer 
to the Local Standard of Rest (LSR) (IReid et al.ll2009l ). Instead, for the astrometric data 
it turns out that only an approximate realization of a common relative reference frame is 
possible. It implies that the exact definition of the coordinates is a matter of each data 
analysis in such a way that simply merging: tw o different sets of astrometric positions would 
yield incorrect results ( IGillessen et al.ll2009bl ). 



In our calculation, we will proceed as follows. For the sake of generality, let us assume 
that an explicit, analytical expression is available for a given observable Y in such a way 
that it is function of all or some Keplerian orbital elements, i.e. Y = Y(f, {k}), where / is 



iSchodel et al.l ( 120021 ) introduced in the scientific literature for the first time such a term 
for designing the pericenter u in the case of a BH. It comes from the Latin word u niqer" , 
meani ng "black". Before, the term "perimelasma" was used by the physicist G.A. lLandis 
(120011 ) in a science-fiction tale of him. Actually, the Greek word "/j,e\a<;" means "black" , 
but, on the other hand, "melasma" (Chloasma faciei) is a medical term denoting a skin 
affection which consists of a tan or dark skin discoloration An earlier and somewhat more 
adequate term is "peribothron" , c oming from the Greek word u [369poq n for "pit" . It is due 
to a suggestion by W.R. Stoeger to lFrank fc Reesl (119761 ). After all, BHs may not be entirely 
"black" since, under certain circumstances, they may radiate, so that a term recalling their 
deep gravitational well appears more suitable to characterize them. Let us note that the 
definition of BHs as "frozen stars" by the former Soviet scholars may have lead to the word 
"peripekton" , from the Greek word "vrr/Kro^" for "frozen" . 
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the true anomaly, and k denotes the ensemble of the Keplerian orbital elements explicitly 
entering Y apart from the mean anomaly M. Then, we straightforwardly compute its 
secular variation as the sum of two parts. The first one is purely Keplerian, and it vanishes 
over one orbital period Pb- The second one is due to the non-Keplerian variations of all 
the orbital elements induced by the dynamical perturbation considered. The total result is, 
thus, 



dY\_ (AY) _ / 1 N 



dY df dM ^dYdK 



dt 

w 



df. (2) 



df dM dt ^ 8k dt 

In it, dM/dt and dn/dt are the instantaneous variational of the Keplerian orbital elements 
computed with, e.g., the Gauss variation equations and evaluated onto the unperturbed 
Keplerian ellipse, while df /dM and dt/df are the usual Keplerian expressions for such 
derivatives: see eq. (ITU1) and eq. ffTTj) below. 

The paper is organized as follows. In Section [2] we deal with some kinds of both classical 
and relativistic perturbing accelerations. The long-term effects caused by them on the 
radial velocity are analytically worked out in Section [3J In Section H] we perform numerical 
calculations by using the S2 star, and confront them with the present-day measurement 
accuracy. Section [5] is devoted to the summarizing our findings. 



2. The perturbing accelerations 

2.1. General treatment and overview 

Here we deal with a generic perturbing acceleration A induced by a given dynamical 
effect which can be considered as small with respect to the main Newtonian monopole 
^Newton = GM,/r 2 , where G is the Newtonian constant of gravitation and r is the relative 
star-SBH distance. The stars orbiting the SBH are assumed test particles: their masses are 
about m* < 10~ 5 M., and relativistic corrections to their internal structures are assumed to 
be too small to yield noticeable effects on their orbital motions. 

First, A has to be projected onto the radial, transverse and normal orthogonal unit 
vectors R, T, N of the co-moving frame of the test particle orbiting the central body 
acting as source of the gravitational field. T heir components, in cartesi an coordinates of a 
reference frame centered in the primary, are flMontenbruck fe Gillllioooh 



cos fl cos u — cos A sin Vt sin u 
R = | sin f2 cos u + cos A cos il sin u ] (3) 
sin A sin u 



2 Actually, dM / dt is the sum of the Keplerian mean motion n and a non-Keplerian term, 
as we will see later. Its Keplerian part yields from eq. ([2]) the Keplerian variation of Y . 
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— sin u cos Q — cos A sin Q cos u 

— sin Q sin u + cos A cos Vt cos u 

sin A cos u 



(4) 



N 



(5) 



sin A sin Q 
— sin A cos Q 
cos A 

In eq. ([3])-eq. (jSJ), fi,w,A are the longitude of the ascending no dfl the argument of 
pericenter, reckoned from the line of the no defl and the inclination of the orbital plane 
to the reference {xy} plane, respectively. In general, we will choose the unit vector p of 
the line-of-sight, pointing from the object to the observer, to be directed along the positive 
z axis, so that the {xy} plane coincides with the usual plane of the sky. With such a 
choice, corresponding to the frame actually used in data reduction fjGhez et al.l 120081 ) . A 
is the familiar i, and Q is an angle in the plane of the sky counted from the mean vernal 
point at J2000 epoch along which the reference x axis is customarily chosen; it is such a 



node w hich is actually determined from the observations fjGhez et al.l 120081 ; iGillessen et al. 



2009aU bh. In other cases, in order to compute more easily certain dynamical perturbations, 
we will orient our frame with the z axis aligned with the central body's proper angular 
momentum L, so that the reference {xy} plane will coincide with the body's equatorial 
plane. In t his case, Q will be an angle lying in such a plane, and it is not the one released in 
literature ( Ghez et al. 2008 : Gillessen et al.l 2009a b). Moreover, u = / + u is the argument 
of latitude. Subsequently, the projected components of A have to be evaluated onto the 
Keplerian ellipse 

p = a(l-e 2 ), (6) 



P 



1 + e cos / 

where p is the semilatus rectum and a, e are the semi-major axis and th e eccentricity, respec- 



tively . The cartesian coordinates of the Keplerian motion in space are ( iMontenbruck &: Gill 
2000h 



x = r (cos cos u — cos A sin Vt sin u) 



r (sin Q cos u + cos A cos Q sin u) 



r sin A sin u, 



(7) 



3 It is an angle in the reference {xy} plane from the reference x direction to the line of 
the nodes. 

4 It is the intersection of the orbital plane with the reference plane {xy}. 
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while the cartesian components of the velocity can be obtained as 



v z 



dxd£ 
df df 

dydt 
df dt' 

dzd£ 
df dt' 



in which df/dt is given by eq. ffTTj) . 

Then, Ar, At, An are to be plunged into the right-hand-s ides of the Gauss equa tions 
for the variations of the Keplerian orbital elements. They are ( IRoyll2005l ; ISoffell 119891 ) 



dt 



de 
dt 

dX 
dt 

<in 

dt 



nVl— e 2 



[Arc sin / + At (-)] , 
{ Ar sin / + At [cos/+ \ (l - r)] } 



— AjvQ sinn, 



n.a\/l— e 
1 



(9) 



da) 
dt 



cosA^ + 



dt nae 



-A R cosf + A T (1 + JJsin/ 



dM 



/i 



nae 



■4cos/ + A T 1 



sin/ 



where n = i/GM/a 3 is the Keplerian mean motion related to the orbital period by 
n = 2n/P h . 

As explained in the Introduction, the right-hand-sides of eq. (Q, computed for the 
perturbing accelerations of the dynamical effect considered, have to be inserted into the 
analytic expression of the time variation dY/dt of the observable Y of interest which, the n 
must be averaged over one orbital revolution according to eq. 02]) by means of (iRoyl 120051 ) 



and 



df 



dt 



Vl-e 2 dM, 



3 2\3/2 



n(l + e cos /) 



;df. 



(10) 



'ir. 
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2.2. The effect of general relativity 



In its slow-motion and weak-field approximation, GTR predicts that a slowly rotating 
central body of mass M and proper angular momentum L induces two kinds of small 
perturbations on the otherwi se Keplerian orbi tal motion of a test particle. The largest one is 
dubbed gravitoelectric (GE) (jMashhoonll2007l ) . and depends only on the mass M of the body 
which acts as source of the gravitational field. It is responsible of the well-known anomalous 
secular precession of the perihelion of Mercury of 43.98 arcsec cty~ x in the field of the S un. 
There is also a smaller perturbation, known as gra yitomagnetic (GM) (jMashhoonl 120071 ). 
which depends on L: it causes the Lense-Thirring ( Lense fc Thirring 19181 ) precessions of 
the node and pericenter of a test particle. In the linearized gravitoelectromagnetic (GEM) 
approximation, the general relati vistic pertu rbing acceleration Aqtr to be added to the 
Newtonian monopole ^Newton is (ISoffellll989l ) 



L GTR 



-E, 



v 



2 - x B 



91 



(12) 



with 



E, 



B r 



GM 



S,\L-3(L 



[( 



AGM 



2 ) r + 4 



v) v\ 



r) r\ 



(13) 



where c denotes the speed of light in vacuum. In eq. (TT2"]) -eq. (JI3"]) E g is the GE field, 
while B g is the G M one . In r egard to B g , in the case of a rotating BH, the existence of 



the horizon in the iKerrl ( 119631 ) metric, whic h describes the spac e time outside it, im plies a 
maximum value for its angular momentum (IBardeen et al.lll972[ iMelia et al.ll200ll ) 



(max) 



MlG 



(14) 



If such a limit is actua lly reached or not by astrophysical BHs depends on their accretion 
history (IBardeenl Il970l ). In fact, recent measurements of the spin of the SBH in Sgr 
A*obtained in the context of discseismology by mea ns of newly detec ted quasi-periodic 
oscillations (QPOs) of radio emission point towards (IKato et al.ll2010l ) 



L.c 

MIG 



0.44 ±0.08. 



(15) 



Genzel et al.l (120031 ) showed that for Sgr A*the spin parameter is 0.52 ± (0.1, 0.08, 0.08) or 
larger, if a QPO observed from Sgr A*in 2003 is of dynamical origin. On the other hand, it 
must noted that X-ray QPOs oscillations are rather disputed. Actually, in the case of the 
S stars orbiting the SBH in Sgr A*, while the weak- field approximation is acceptable since 



for S2 



(U) GM. 1 1 



GM. 



4 x 10" 



(16) 
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the slow-motion approximation may be, in principle, less adequate. Indeed, the speed of, 
say, S2 at perinigricon is as large as 2.6% of the speed of light, while at aponigricon it is 
0.1% of c. Thus, in the dynamical equation of motion of ea. f[T2|) higher order relativistic 



corre ctions should be, in principle, taken into account (IMashhoonl l2005t ICapozziello et al. 



20091 ). Anyway, as we will see, they induce negligible consequences on the dynamical effects 



we are interested in, given the present-day level of accuracy in S stars spectroscopy. 



2.2.1. The gravitoelectric, Schwarzschild-like perturbation 



By defining 



GM 



(17) 



the R — T — N components of the general relativistic GE perturbing acceleration E g are 
f lSoffellll989[ ) 



.4 



(GE) 



n 2 Tlg 
(1-e 2 ) 3 



R 
/l(G E ) 



1 + e cos /) 2 (3 + 2e cos / - e 2 + 4e 2 sin 2 /), 



0fy (1 + e cos /) 2 4e sin f(l + e cos /), 



(18) 



.4 



(GE) 
N 



0. 



Note that [7Z g ] = L, so that [n 2 lZg] = L T 2 . Moreover, eq. fTTSl does not depend on the 
inclination of the orbit to the plane of the sky. 



2.2.2. The gravitomagnetic, Lense-Thirring-like perturbation 



The R — T — N components of the Lorentz-like general relativistic GM perturbing 
acce leration ind uced by the rotation of the central body with proper angular momentum L 
are flSoffellll989h 



with 



.4 
.4 



H 

(GM) 
T 

(GM) 
N 



rjg cos \P(1 + e cos /), 
— r] g e cos ^ sin /, 

Vg sin vD(l + e cos /) [2 sin u + e (f^ff 



(19) 



Vg 



a 2(! _ e 2)7/2 



1 + e cos fy 



(20) 
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and 

^ = (2D 



c 2 



Note that [Q = L 3 T -1 , so that [%] = L T~ 2 . In eq. (USD * is the angle between the 
orbital plane and the equatorial plane of the central body, i.e. L has been assumed directed 
along the positive z axis so that the reference {xy} plane coincides with the equatorial 
plane of the central body (A — > \&). For equatorial orbits, i.e. for \I/ = 0, A^ M ^ = and 
A^ M ' ^ 0, A^ M ^ ^ 0. Instead, for polar orbits, i.e. for \1/ = 90 deg, only the normal 
component does not vanish. 



2.3. The quadrupole mass moment 



The external gravitational field of a rotating body undergoes departures from spherical 
symmetry because of the distortion of its shape due to the resulting centrifugal force. 
An oblate body of equatorial radius R e and adimensional quadrupole mass moment J 2 
affects the orbital motion of a test particle with a non-central perturbing acceleration 
f lCunningham]|l970l : IVrbildl2005h 



3J 2 R 2 e GM 
2V 4 



5(f-L) 



r + 2(r- L)L 



(22) 



where L is the unit vector of the body's angular momentum, directed here along the 
positive z axis so that the equatorial plane is the reference {xy} plane (again, A — > ty). 
According to eq. ©-eq. (JSJ) and eq. (J7J), the R — T — N components of eq. fl22|) are 



WO 



.4 



AJ2) 

l\.rp 



A 



(J2) 

N 



A( J 2) . J< - 

A {-h) . 



3n 2 R-J, 
'8a(l-e 2 ) 



t (1 + e cos ff (1 + 3 cos 2^ + 6 sin 2 * cos 2u) 



n R e .h_ (! + e cos ff sin 2 y gin 2u ^ 

1 + e cos f) A sin 2\1/ sin u. 



2a(l-e 2 ) 

3n 2 R 2 e J 2 
2a(l-e 2 ) 4 



(23) 



Note that [n 2 R 2 aT l ] = L T -2 . For ^ = 0, i.e. for equatorial orbits, only the radial 
component is not zero. For polar orbits, i.e. for $ = 90 deg, the normal component 
vanishes, contrary to the radial and transverse ones. 

Also a rotating BH should have a quadrupole mass moment, so that it affects the 
orbital motion of a distant test particle with a perturbing acceleration analogous to that 
of eq. ( l22l) . It is customarily to i ntroduce a dimensional quadrupole parameter Q 2 , 



[Q 2 ] = L 5 T~ 2 , in such a way that ( 1Willll2008l ) 



J 2 R 2 e GM -> Q 2 



(24) 
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throu ghout eg. (1221) and eg. (|23|). A ccording to the "no-hair" or uniqueness theorems of 
GTR (jChrusciell Il994j ; lHeuslerlll998l ). an electrically neutral BH is completely characterized 
by its mass M. and angular momentum L, only. As a consequence, all the multipole 
moments of its external spacetime are functions of M. and L,. In particular, the quadrupole 
mass moment is 



Ql 



LjG 
c 2 M, 



-xl 



G 3 M? 



(25) 



Thus, in the case of the SBH in GC eq. ( fTBl yields 



\Q'\ 



3.585 x 10 45 m 5 s" 2 . 



(26) 



In terms of the adimensional coefficient J2, by assuming the Schwarzschild radius r g for the 
equatorial radius R e of the BH, eq. ( 1261) would correspond to 



\J'\ = 4 x 10" 



(27) 



Just for the sake of a comparison, for the Sun, Jupiter and the Earth we have = 2 x 10 



from helios eismology (IMecheri et al.l 120041 ) , J, 



r (Jup) 



spacecraft (IJacobsonl 120031 ) . J® 
( lBruinsmall2010l ). respectively. 



1.46 x 10 from the flybys of some 
1.08 x 10~ 3 from the dedicated GRACE spacecraft 



2.4. Inner diffuse mass distribution 



It is well recognized that, in addition to the dynamical effects considered so far directly 
related to the SBH, also the impact of a diffuse cluster of non-luminous ordinary mattecl 



around the BH due to massive remnants of various kinds (Morris! Il993l ) should be taken 
into account. Indeed, reasoning in terms of the perinigricon, such an extended material 
component would induce a retrograde precession cu, which may overwhel m the general 
relativistic GE one for certain values of its mass f lRubilar fc Eckartll200l[ ). Let us, now, 
work out in detail its dynamical effects. 



Following iRubilar & Eckartl f l200llUMouawad et al.l ((20051), iGillessen et al.l fl2009bh . we 
adopt a Plummer density profile 



Qdm{r) 



3/iM. 

Audi 



r 2\ "5/2 



(28) 



5 In principle, one should also take into account t he issue of non-baryonic dark matter 
cusps in relation to the merging history of galaxies (IMerritt fc Milosavljeviq 120021 ). It is 
beyond the scopes of this work. 
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where the core radius is (IGillessen et al.ll2009bl ) 

d c = 15 mpc, 



(29) 



in agreement with the observed light profile (IMouawad et al.ll2005l ). while n, the mass 
parameter, is the ratio of the total extended mass M T at a given distance r to the central 
point mass. IGillessen et al.l (l2009bl ). with a fit involving S2 able to probe the mass enclosed 
between its aponigricon and perinigricon, yield 



/i < 0.04. 



(30) 



It should be mentioned that IMouawad et al.l (120051 ). with the position of Sgr A*as a free 
discrete input parameter, provided for the first time an upper limit of 



fx < 0.05. 



From the Poisson equation 



1 d 



,8U 



dr 



AttGq 



(31) 



(32) 



for the gravitational potential U, written in spherical coordinates since U = U{r) in view of 
the fact that g = g(r), the perturbing acceleration 



easily follows. It turns out 



L(dm) 



^(dm) 
R 



dU 



(dm) 



GM.fi 



(dm) . 

r 



Or 



(33) 



717 



-3/2 



^(dm) 
.4 



T 
(dm) 



0. 
0. 



(34) 



3. The radial velocity 

The basic observable in spectroscopic studies of binary systems is the radial velocity 
v p , i.e. the component of the velocity vector v of one of (or both) the system's partners 
along the line-of-sight whose unit vector p has been assumed directed along the z axis. 
Concerning the S stars orbiting the Galactic SBH, let us mention the possibility that some 
of them may actually be binaries themselves. In this case, additional "noise" to the radial 
velocity data would be introduced. Its expression for unperturbed, Keplerian elliptic 



13 



orbits, up to the velocity of the system's center of mass Vo, can be obtained by using the 
z compone nts of eq. ([3]) and eq. ((!]) with A — > i, and recalling that, for a Keplerian orbit 
(fik^l2()().',h. 



v = vrR + vjT 



na 



esin fR+ (1 + ecos/)T 



The result is 



■v. 



K [e cos to + cos(/ + u)] , 
where 2K is the total observed range of radial velocity defined by 

na sin i 



K 



(35) 
(36) 

(37) 



Note that eq. (136]) an d eq. (1371) agree with the expressions given by, e.g., iBattenl ( 120011 ) and 



Padmanabhanl ( 120101 ). 



Perturbing dynamical effects affect the radial velocity as well by inducing, in principle, 
a non-vanishing net radial acceleration over one orbital period. It can straightforwardly be 
worked out from eq. ([2]) with Y — > v p by noting that, in this case, the perturbations of all 
the six Keplerian orbital elements are involved. In this respect, a special care is required for 
the parameter i entering eq. ( 1371) . It is the angle between the unit vector i of the orbital 
angular momentum and the unit v ector p of the line-of-sight pointing towards the observer. 
From the spherical law of cosines (IGellert et al.l Il989l ; IZwillingerl 119951 ) 



cos B = sin C sin A cos b — cos C cos A 
with the identification^] A — > B — > n — i, C — > I, b — > n — f2, it turns out 

cos i = sin \1/ sin / cos Q + cos \1/ cos /, 



(3f 



(39) 



where \1/ is the angle between £ and the unit vector L of the central body's proper angular 
momentum, I is the angle between L and p, and Q is the longitude of the ascending node 
defined from 

sin \1> sin / cos — (L X 1) • (L X p). (40) 

From eq. (|39|) it turns out 



sini 



It) 



fsin \I/ cos / — cos \I/ sin / cos fi) h sin ^ sin / sinf2 ( — ) 

k 1 dt \dt J 



(41) 



As eq. (EJ) shows for A — > ^, only those perturbing accelerations, like the general relativistic 
GE one, having no out-of-plane components An leave i unaffected since both ^ and Q may 



6 In such a way, results to be prograde with respect to the orbital motion, i.e. Q follows 
it, coherently with the definition of the longitude of ascending node. Moreover, it lies in the 
equatorial plane. 
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change due to A^. Thus, in the following expressions for the effects due to the general 
relativistic GM field and Q2 it must be recalled that Q is intended to be referred to the 
equatorial plane, not to the plane of the sky. It must also be noted that both / and \I/ are 
unknown. 



3.1. The effect of general relativity 

General relativity dynamically affects the radial velocity of non-circular orbits by 
causing averaged long-term variations of it. Indeed, an exact calculation in e with eq. ffT8j) 
yields 

/ . rrpn / 9^ \ 15e(l + e 2 ) sinz sin a; . „ , 

{vf E )> = (nX) 3 (1 J e2)5/2 > (42) 

for the Schwarzschild-type, gravitoelectric component, while the secular acceleration due to 
the Lense-Thirring-type, gravitomagnetic terms of eq. (1T91) i^j] 

W GM) > = (^) 4(1 ^ e2)2 [Ve(i, /, 0, *) cos a; + V.(i, J, 0, ) sin w ] , (43) 

with 

V c = 11 cot i sin / sin sin Q, 

y s = {cos n sin 2/ (sin * - sin 3*) - (44) 

- sin 2 * cos * [cos 21 (3 + cos 2Q) + 2 sin 2 fi] +104 sin 2 i cos * } . 
Let us recall again that / is the angle-unknown-between L and the line-of-sight , and Q lies 



in the equatorial plane of th e source: for it the values determined in literature (IGhez et al. 
20081 : iGillessen et alJ E?>n9al fbl) cannot be used because they refer to the plane of the sk}o 



Also is unknown. Anyway, in order to give an order-of-magnitude evaluation of the 
Lense-Thirring long-term radial acceleration we will simply look at the multiplicative 
dimensional factor in front of the square brackets in eq. (I43p . For e — > both eq. fj4*2l) and 
eq. (I43p vanish. The general relativistic effects are non-vanishing either for edge-on (i = 90 
deg) or equatorial orbits (\& = 0), with eq. (|42p which is independent of contrary to eq. 



7 It is assumed that the orbital motion is prograde with respect to the BH's spin; a — 
sign would occur in the opposite case. 

8 In principle, spherical trigonometric formulas relate both the nodes Q e and Q s , where e 
and s denotes the equator and the sky. However, the knowledge of the angle, in the plane of 
the sky, between the equatorial plane and the reference {xz} plane would be required. 
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Some preliminary and approximate calculation of the effects of both the general 
relativist ic GEM dynamical effec ts on the velocity of S stars orbiting the SBH in GC can be 



found in iKannan fc Saha I (120091 ) ; they do neithe r deal with the d irectl y measurable radial 
velocity nor with its possible variations. See also lAngelil &: Sahal (120101 ) for GEM effects on 
the traveling photons paths. 

Special r elativistic effe c ts on the radial velocity related to the Doppler effect have been 
considered by IZucker et al.l (120061 ) ; anyway, they do not involve net variations of v p over one 
orbital revolution. 



with 



3.2. The quadrupole mass moment 

The BH's oblateness causes an averaged long-term variation of the radial velocity only 
if the orbit is elliptic. Indeed, it turns out 

k {j2) ) = - m) — — %m — ^ j ' ^ *) c ° suj + j ^ j > ^ n > *) ^ > 

' V a J 32 (1 — e 2 ) ' sin % 

(45) 

J c = 10(1 — e 2 ) cosi sin /sin 2\I> sin f2, 
J s = 2 (1 — e 2 ) cos % sin 2^ (cos /sin \l/ — sin / cos \l/ cosO) + 

(46) 

+ sin 2 i [7 + 47 cos 2* + sin 2 ^ cos 2a;- 

- J (259 + 429 cos 2^ - 44 sin 2 * cos 2w) 

It is an exact result in e, and vanishes in the limit e — > 0. Note that, according to eq. ( 146]) . 
eq. (1451) vanishes neither for equatorial orbits (\l/ = 0) nor for edge-on configurations (i = 90 
deg). Also in this case, Q refers to the equatorial plane. 



3.3. The diffuse inner dark matter 



The averaged long-term effect of the diffuse inner dark matted on v p can be worked 
from eq. ( )34|) . An approximate calculation with 



r 2 V 3/2 3r 2 



9 The effects of non-baryonic dark matter on the perinigricon precession have been con- 
sidered by, e.g., IZakharov et al.l ( 120071 ) . 
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and 
yields 



(1 + ecos/) 4 w 1 - 4ecos/, (1 + ecos/) 5 ^1 — 5ecos/, 
GM.a\ e(l -e 2 ) 3 / 2 /isinz 



dl 



32 - 5e J (2 + 3e 2 ) + 6(1 - ey(5e 2 - 11 



(48) 



sinw. 

(49) 



Note that eq. (149]) vanishes for e — >■ and for face-on (i = 0) orbital configurations. 
Concerning the validity of the approximation of eq. (147|) . it actually holds for S2. Indeed, 
its orbital parameters and eq. ( 129]) for d c tell us that 



0.002 < — < 0.4. 

di 



(50) 



4. Numerical evaluations and confrontation with the measurement accuracy 

By using the know n orbital parameters of S2 along with the as sociated uncertainties 
(ICillessen et al.ll2009bl ). the known mass M. of the SBH in the GC dOillessen et al.ll2009bl ) 
and the values for its angular momentum and quadrupole mass moment from eq. ( fT5]) and 
eq. ([26]), it turns out tha^ 



4 GE) ) 



(dm) 



(GM) 



5.2 ±2.6) x 10~ 5 m s" 



3.8 x 10" 6 m s~ 2 



oc 1.3 x 10" 8 m s" 



(51) 



vf 2) ) oc 1 x 10- 1 



m s 



Concerning a possible measurement of a net change in the radial velocity of S2 after it 
completed one full orbital revolution, no empirical evidences for it exist to date in literature, 
at least to the knowledge of this author. Anyway, measurements of the magnitude of the 
three-dimensional acceleratio n of S2 after 2 y ears (1997-1999) are available; its accuracy 
amounts to 4 x 10" 4 m s" 2 flGhez et all 120001 ) . By_ assuming an unc ertainty of about 15 
km s" 1 in measuring the radial velocity of S2 ( Gillessen et al.l 2009a ). an overall accuracy 



10 



The figures for (i^ M ^ and (i^ 2 ^ refer to the dimensional multiplicative factors in 



front of the square brackets in eq. ( 143]) and in eq. (|45 
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of the order of 2.4 x 10 m s in (v p ) may be assumed in future over an observational 
time span At = 20 yr. Actually, it must be pointed out that the currently available radial 
velocity measurements do not yet cover one full orbit revolution for S2. Indeed, the first 
radial velocity data points are from 2000, then 2002; they are more densely sampled from 
2003 onwards (S. Gillessen, private communication, August 2010). 

It has to be pointed out that the total accuracy reachable in the changes in the radial 
velocities is actua lly impacted by the uncertainty in LSR itself. Indeed, as explained by 
(IGhez et al.ll2008l ). to obtain the radial velocities with respect to the LSR, each observed 
radial velocity has to be corrected for the Earth's rotation, its motion around the Sun, and 
the Sun's peculiar motion with respect to the LSR (nominal value U & = 10 km s _1 , directed 
radially inwards, ( Dehnen fc Binney 1998h ) . Since the LSR is defined as the velocity of an 
object in circular orbit at the radius of the Sun, the Sun's peculiar motion with respect to 
the average velocity of stars in its vicinity should give the Sun's motion toward the GC. 
In all such a mac hinery, also the rotation speed Oo of LSR enters; recent evaluations by 
Reid et all (120091 ) yield an uncertainty of the order of 16 km s _1 corresponding to a future 



uncertainty of 2 x 10 5 m s 2 over At = 20 yr. Moreover, also the motion of the SBH itself 

In pa rticular, the uncertainty in its radial 



should be taken into account (IGhez et all 12008 

-l 



(jGouldl 120041 ) implying a limit in the accuracy in 
20 yr. Such limiting factors should be taken into 



velocity can be evaluated to be 2 km s" 
(v p ) of about 1 x 10~ 6 m s~ 2 over At = 
account when future improvements in measuring rad ial velocities are discu ssed. Searches 
for pulsars orbiting the SBH are currently underway (jMacquart et al.ll2010l ); their discovery 
may yield orbiting probes with a better accuracy in their radial velocity changes. 

These considerations show that while the gravitomagnetic and the quadrupole effects 
are far from being directly detectable in such a way, the gravitoelectric trend lies just at the 
edge of the measurability capabilities. The effect of the diffuse dark matter inside the orbit 
of S2 is one order of magnitude smaller than the general relativistic GE one. Incidentally, 
such figures indicate that higher-order corrections to the computed effects of eq. (I42l) -eq. 
(14"3"|) in the linear GEM approximation due to the relativistic motion of S2 can be neglected 
at the moment. 

Finally, in view of the possible discovery in the near future of stars closer than S2 we 
plot in Figured] the magnitudes of the effects considered as functions of the semi-major axis 
for different values of the eccentricity. 



5. Conclusions 

One of the directly measurable quantities of the system constituted by the S stars 
orbiting the Supermassive Black Hole located at the center of the Milky Way in the radio 
source Sgr A*is the radial velocity. Given that S2 has already completed one full orbital 
revolution, with a period of 15.9 yr, since its discovery, and in view of the possible detection 
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Fig. 1. — Maximum values of the long-term time variations (v p ), in m s -2 , as a function of a 
(500 au < a < 1031.69 au) for different values of the eccentricity: e = 0.1 (red dash-dotted 
line), e = 0.3 (blue dashed line), e = 0.6 (light blue dotted line), e = 0.9 (yellow continuous 
line). For the the spin L, and the quadrupole Q\ of the SBH we used the values of eq. ([TBI) 
and eq. f[26l) . respectively. 
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in the near future of other stars and pulsars with shorter orbital periods, we looked at the 
cumulative, long term time variations of the radial velocity caused by several Newtonian 
and Einsteinian dynamical effects. They are both the general relativistic Schwarzschild 
and Kerr-like components of the spacetime metric, the quadrupole mass moment and the 
diffuse dark mass distribution made by stellar remnants enclosed within the star's orbit. 
We analytically worked out the long-term variations in the radial velocity induced by them 
by finding non-zero effects for all of them. We used S2 for numerically computing their 
magnitudes. They are 8 x 10~ 5 m s~ 2 (Schwarzschild), 4 x 10 -6 m s~ 2 (dark matter), 
1 x 10 -8 m s -2 (Kerr), 1 x 10~ 10 m s~ 2 (quadrupole), respectively. In computing the 
general relativistic variations of the radial velocity, we remained within the post-Newtonian 
regime by neglecting relativistic corrections of higher order in the equations of motion. The 
figures for the the Kerr and quadrupole effects have been computed by using the latest 
determinations of the angular momentum parameter of the Galactic black hole, and in the 
"no-hair" hypothesis, respectively. For the dark matter distribution we used a Plummer-like 
mass density profile. By assuming a present-day uncertainty of about 15 km s _1 in the 
radial velocity measurements, its time changes may be detected in the future at a 10~ 5 
m s~ 2 level over an observational time span of 20 yr; at present, radial velocity data cover 
just 7 yr. Even if such evaluations will turn out to be not too optimistic, a detection of 
the Kerr and the quadrupole-induced cumulative changes of the radial velocity seems to be 
unfeasible. 
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